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PROBLEM OP GAS ACCRETION ON A GRAVITATIONAL CENTER 

V. A. Ladygin 


Examined in the study is a method of approximated solution 
of the problem of gas accretion on a rapidly moving gravitational 
center. The solution is obtained in some vicinity of the axis 
of symmetry in the region of potential flow. Calculations on 
a computer showed the effectiveness of the given method. 

Introduction 

The solution of the problem of stationary gas accretion on 
a moving gravitational center simulates the movement of a 
substance in interstellar space in the vicinity of a black hole. 

A detailed picture of gas accretion on a black hole is of 
interest in connection with the problem of observation of black 
holes. 

The qualitative study of such accretion, as well as two- 
dimensional numerical calculations of this problem, are available 
in studies Ll] - 03* The study of self-modeling solutions, which 
may represent asymptotics of the flow of gas near a gravitational 
center, was carried out in study L5 U • 

In the present study, the system of equations of two- 
dimensional gas dynamics, which describes gas accretion, in 
contrast to studies L .23 - 03 » is solved in an approximate 
manner, by means of expansion into a linear series, according to 
one of the independent variables (angle 0), and by means of 
" abridging” of the obtained infinite system of common differ- 
ential equations. 


♦Numbers in the margin indicate pagination in the foreign text. 



Although only the region of potential flow of the gas is 
studied in the present study, this method may he used also for 
a nonpotential flow behind a shock wave. 

1. Formulation of the Problem and Basio Equations 

Studied herein is the steady-state axisymmetrical flow 
of an ideal polytropio gas, devoid of viscosity and thermal 
conductivity in a gravitational field of a material point of 
mass M. 

At infinity, the approach stream is assumed to be homo- 
geneous and supersonic. 

p-p^const t p-*R.«const t 

, /2.r/ 

V^CUT15t>0, * l ~*" 


With these boundary conditions, the flow prior to the 
shock wave is potential and isentropic, and, in a spherical 
system of coordinates, it is described by the system of 
equations t , . k 




-energy integral, 


/S2.3/ 


-condition of potentiality, 

4»^(pt’V,) + . Y1! |i rs -^(pV4sinS)-o /i . «/ 

-equation of continuity, 

P-(r -w* /2>6/ 

-equation of state, 


where V< , V# ,j0 , and P are the radial and angular com- 
ponents of velocity, the density, and the pressure of the gas, 

G is the gravitational constant, y is the indicator of the 
polytropic curve . 


The constants K and Y are determined with the boundary 



conditions (2.1). 

K-fF&K . '*'* ^ ' (Vntc ■ nA/ 

Through transformation of the analog 

7 , p^fLpi, Vt - CVi. , V fc (J4, P I'flft , 

where r - %>»• is the speed of sound at infinity, the five- 
parameter problem (2.1)-(2.6) (parametersi y»G*M, P® , V®») 
is reduced to a two-parameter problem (parametersi y and M*« 

is the Mach speed at infinity) . 

C*» 


Therefore, without bounding the generality, one can assume 
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42. V 


On the strength of (2.4), the expression 

-ipsinOl/, +rps tnOV t dJ 

is a complete differential of some function S(1 ,6) of the 
current lines, and, consequently, 

d|$ «*-T$ift0Vi *| 

dff** x'stnOVt ! 

/2.2/ — Ml. ; 

/&*s/ — — ** v#) - » o ’ 


n >if 

A*..o; 

A.il’ 


•j 

where d=jp is the specific volume. 


2. Description of the Method of Solution of the Problem 

We will derive the formulas for the approximated solution 



of the system (2.8)-(2.11) . Since the flow is axisymmetrical, 
then 




W© will assume that the functions d, V* , V* , S are 
analytic according to 0 in some vicinity of 9*0. Then, d, 

VW , S are ex* ' 'ided into an exponential series according to 
even powers of 0 , and V* , according to uneven powers of 6 , i.e. « 


dMi-CjUv 0 “, r . 

S(t.w.£s„tt> e u . v.k.sj-I!, p„.«> a 


W 


One can think that S* ('t)*0, since the function S(t,8) is 
determined with an accuracy up to the additive constant, and is /6 
constant along the trajectory, while the axis of symmetry 0=0 
is the trajectory of the particles. 

We will substitute the expansion (3*2) into the system 
(2.8) -(2.11) , and group the terms with identical powers of 0. 

We will obtain the system of relationships » 


/i *ti,- 

/:\9/ 

A’.IC/ 

/a. ii/ 


* f 


/k*1 .2, 


■ it/ 1 


/a. 3/ 
/3.4/ 

/a.e/ 

/3.6/ 

/3.?/ 


where a* , , h t ' are the coefficients of the expansion into 

an exponential series, respectively, of the functions sin 0, 
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v’+v* , 




4m0 0 u , - iSt- tyr ' 

/ K * 0 1 1 1 % * */ 1 

d' r ~h*£ih< r 


/3.6/ 

/B.9/ 


We will expand the equality 

dh t frit d* 

into an exponential series according to 0 . 


We will obtain a system of relationships between the 
coefficients Ify} and id^}. 


h* y . 

} V /•, u' ~*vl ) / l t h 


/3.I0/ 


We will transform the system (3*3) -(3.10) . In place of the 
coefficients S aK , we will examine the functions 

/3»n/ 

r a i** h l* / 

* * * f 

which have finite limits with t”** 00 . Using (3*5) and (3.6), 
we preclude the coefficients hi in (3.10). We will introduce 
the new independent variable 


U-fiil 


/3.t2/ 


and also the functions 

* dm / d* 

/k*I*2 , « 4 */ • 

We will obtain an infinite system of common differential 
equations relative to the functions 


^-2 

-r„-F,A2-& 




Z2 
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The functions [T^ k }^ , Iq.*, 3*,, » Q> 8 (square of the speed of 
sound on the axis of symmetry 6*0), which are part of the right- 
hand portion of system (3*13) » as well as the functions 
Ujoc » are sequentially determined through the relationships! 

V,*4fi\ X»~Pt 

Then, if the coefficients \i d , V^.q^ (£<K) and T^ (j<K+l) 
hav' already been calculated, then p^, q^ are determined 
from the system of three linear equations! 



/rue/ 

/h.i?/ 

/3.IP/ 


and the function T,^,,) is calculated according to the 
formula 


z 


H*M i 




C p. ^ 
.*)«{*•! 


/;i ! ’•' 


From (2.1), (2.9), and (3-2), we will obtain the boundary 
conditions for system (3*13) 


with tt -*+*» 


(- i '*' 1 m 

T7k"« v 



/jjv 


/*#*,*, ••/ 


If the system (3*13) -(3*20) is solved, then the coeffici- 
ents d 4K and are determined by the formulas 

d t * cf* , dit, 9 * d®/ 4 ** , /3.2I/ 

/k»I ,2 ‘ 


Z8 
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The circulation® aroording to formula® ( 3 • 13) -C 3 • 19) « 
with small values of X » require® large outlays of computer 
time, since tha right-hand portions of tha system (3*13) are 
unbounded witht^O* Therefore, in the region of u<0, it is 
advisable to make a substitution of the variable s« 


pt < m e H pu-< 

q: 

c: -tf'n . . 


F. r-e^K. , 

, /9,22/ 

vu*«e k u* 


In this case, only equations (3*13) and (3*14) change, 
taking on the form* 


'4*5 


u_t 


i- * 


df’ 

Tm 


in » 
« 


r» -h r i f s 

'vjr t'** t 




r,-> 


/3#83/ 

/3«S<t/ 


For approximated calculation of the firat N coefficients 
of expansion (3*2), we will examine .system (3*13) for k*X,2, . . .N. 
Formulas (3 •!**■) -(3*19) make it possible, in the right-hand portion 
of (3*13) » to express q a ' { p^, V J( * (t<N) and T^* (,/< N+l) in the 
form of functions fromu,/ 3 * 5 , , . . . F A , , . . ,iy For a similar 
recording of the coefficients q„ and V M , «rs will make use of 
formulas (3*16) and (3*17) for h s N, withfV*Q being assumed in 
(3,17). This additional assumption does not contradict the 
boundary conditions 

M*' - »> with u •* ♦** 


3, Discussion of Results 

The obtained closed system &N of common differential 

equations relative to the functions p, , . . . ,p Jt) ^ , Fj , . . . ,F^ 



was integrated numerically on a computer by the Runge-Kutt 
method in the interval -ktniO<u<iOtniO (j.e, . 10~*<7<L0 /£> ) 
for the indicator of the polytropic curve y*5/ 3 and the Mach 
speed M#*2.4. In this case, the boundary conditions with^~*+ w 
were transposed to the point T *10 w Bondi radii* 


The calculation was carried out for N*l, 2, ... ,10,15,20. 

For N>2, within the limits of accuracy of integration 
£=10"* of the system (3*13) » the values of the density (X ) , 
Mach speed M(t) , and modulus of speed V(t) on the axis of 
symmetry 0*0 practically do not depend on the selection of N, 
and are represented in the table. The density, speed, and 
Mach speed approach infinity monotonously withl-Q. 

The solution of system (3*13) may be continued into the 
randomly small vicinity of the pointy *Q, which indicates the 
absence of a departed shock wave in front of the gravitational 
center for a gas with y s 5/3* The density J 3 (t ,0) , modulus of 
speed V(T,0), and function of the current lines S(t,0) were 
calculated approximately according to the formulas 




/4,I/ 

V«(V»+V,-8+ , 


/4J2/ 

S -S,0‘* ..*S u 9‘“ 

/4.3/ 


For N-10, the level lines V, and S are given in figures 
1, 2 V and 3 , respectively. Here, the x-axis is the axis of 
symmetry. The gravitational center is located at the origin of 
the coordinates. Plotted along the axes of the coordinates 
are the distances in Bondi radii (R b sG*M/C^ ) . The gas flows 
into the center from the right. 

The calculation of the level lines ofJ>, V, S for N=5, 10, 

15, and 20 shows that, in the range O<0 <tt/ 2, the values of J> , V, 

and S practically do not depend on N, i.e. . there occurs /lQ 
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convergence of the approximated solution to the potential flow. 

For tr/2 <0<ir* and especially olose to the axis of symmetry 
0*Tf r the approximated values o£j>, V, and a depend strongly on 
N. In this region# there is no convergence to the point 
solution, since the flow is nonpotential. 
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